Optical follow-up of GRB 970508 
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ABSTRACT 

We report on the results of optical follow-up observations of the counterpart of 
GRB 970508, starting 7 hours after the event. Multi-color U, B, V, R c and I c band 
observations were obtained during the first three consecutive nights. The counterpart 
was monitored regularly in R c until ~ 4 months after the burst. The light curve after 
the maximum follows a decline that can be fitted with a power law with exponent a. 
= —1.141 db 0.014. Deviations from a smooth power law decay are moderate (r.m.s. = 
0.15 magnitude). We find no flattening of the light curve at late times. The optical 
afterglow fluence is a significant fraction, ~ 5%, of the GRB fluence. The optical 
energy distribution can be well represented by a power law, the slope of which changed 
at the time of the maximum (the spectrum became redder). 

Subject headings: gamma rays: bursts — gamma rays: individual (GRB 970508) 

1. Introduction 

With the Wide Field Cameras (WFCs; Jager et al. 1995) onboard BeppoSAX (Piro et al. 
1995) it became possible for the first time to rapidly determine accurate (3' radius) GRB positions 
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(Costa et al. 1997a,b; In 't Zand et al. 1997; Heise et al. 1997a); such positions, as it turned out, 
are the necessary ingredient for succesful counterpart searches (see also Takeshima et al. 1997 and 
Smith et al. 1997). The detection of a fading optical counterpart to GRB 970228 provided the 
first arcsecond localization of a GRB (Groot et al. 1997a; Van Paradijs et al. 1997; Galama et al. 
1997a,b). Its coincidence with a faint extended object (Groot et al. 1997b; Metzger et al. 1997a; 
Sahu et al. 1997a), which is probably the host galaxy, suggests a cosmological distance for the 
GRB. A second optical transient (OT) was found for GRB 970508 (Bond 1997). The measurement 
of a lower bound to the redshift {z > 0.835; Metzger et al. 1997b) of this optical counterpart has 
unambiguously settled the cosmological distance scale to GRB sources. 

On May 8.904 UT the Gamma-Ray Burst Monitor (GRBM) on BeppoSAX recorded a 
moderately bright GRB (Costa et al. 1997c), which was also recorded with its Wide Field 
Cameras. The fluences, recorded with the GRBM (40-700 keV) and WFC (2-26 keV), were (1.8 
± 0.3) xl0~ 6 ergcm~ 2 and (0.7 ± 0.1) xl0~ 6 ergcm~ 2 , respectively (Piro et al. 1997a). The 
burst was also recorded (Kouveliotou et al. 1997) with BATSE (Fishman et al. 1989) with a 
total 20-1000 keV fluence of (3.1 ± 0.2) x 10~ 6 erg cm" 2 , and a peak flux density (50-300 keV) of 
(1.66 ± 0.06) x 10~ 7 ergcm~ 2 s -1 . From optical observations of the WFC error box (Heise et al. 
1997b), made on May 9 and 10, Bond (1997) found a variable object at RA = 06 h 53 m 49!2, Dec 
= +79°16 / 19"(J2000), which showed an increase by ~1 mag in the V band. BeppoSAX Narrow 
Field Instrument observations revealed an X-ray transient (Piro et al. 1997b) whose position is 
consistent with that of the optical variable. Extended emission is not associated with this optical 
counterpart (Fruchter et al. 1997, Sahu et al. 1997b, Pian et al. 1997), and Natarajan et al. 
(1997) showed that either GRB 970508 originated from an intrinsically very faint dwarf galaxy or 
that it occured at a large distance from a host galaxy (> 25h^Q kpc). 

We here report on the results of optical photometry of GRB 970508, made between 0.3 and 
110 days after the burst. 

2. Optical Photometry 

We have used the 4.2m William Herschel Telescope (WHT) at La Palma to make optical 
observations in U, B, V, R c and I c of the WFC error box both with the Prime Focus Camera (PF, 
f.o.v. 9^0x9^0), and the Auxiliary Port Camera (AUX, f.o.v. 0'.9 radius). All WHT PF exposures 
lasted 600 seconds, except the one made on Aug 26.9 UT which lasted 1800 seconds. The WHT 
AUX exposures were 3600 seconds each. As part of the Kitt Peak National Observatory Queue 
observations, observations were also made with the 3.5m WIYN telescope, using the WIYN Imager 
S2KB at the Nasmyth focus. Ten 6f8x6^8 R-band images were obtained between May 9.15 and 
9.23 UT at several pointings to cover the whole error box. All exposures lasted 300 seconds. A log 
of the observations is shown in Tab. [p. 

We have corrected the U and B band images for the non-linearity at low counts of the 
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L0RAL2 CCD Chip[5 used for the May 10, 11, and 12 WHT observations (this effect is negligible 
in the V, R c and I c observations). We made a photometric calibration (U, V, R c and I c ) using 
WHT observations on May 9.99 UT of the Landolt (1992) fields PG1047+003 and PG1530+57. 
For V, R c and I c we find good agreement with the results of Sokolov et al. (1997). To include a 
B band calibration, and avoid small calibration differences with respect to Sokolov et al. (1997) 
we used the B, V, R c and I c reference stars (1,2 and 4 in Fig ||) as given by Sokolov et al. (1997) 
(their star 3 was close to saturation and therefore not used). Stars 2, 5, 6 and 7 were used as 
secondary photometric standards for the U band (see Fig ||). We have corrected for a U band 
atmospheric extinction of 1.3 times the nominal extinction (as Carlsberg Automatic Meridian 
Circle V band extinction measurements indicate) 12 . 



3. The R c Band Differential Optical Light curve 

In Fig. [I] we present the R c band optical light curve, where we have included results of 
differential photometry (relative to local reference stars; filled symbols) and calibrated photometry 
(open symbols). The differential magnitudes in Tab. |l[ were determined with respect to stars 1, 2 
and 4. The other differential data are from Sokolov et al. (1997; determined with respect to star 
1, 2, 3 and 4), Pedersen et al. (1997), Garcia et al. (1997), Schaefer et al. (1997) and Chevalier 
et al. (1997). The magnitudes in the last four data sets were determined with respect to star 4 
only. In the absolute calibration of these differential data sets we used the R c magnitudes of these 
reference stars as given by Sokolov et al. (1997). Note that this involves a correction to the data 
of Pedersen et al. (1997), Schaefer et al. (1997) and Chevalier et al. (1997) by -0.21 magnitudes 
for the difference by that amount in R c for reference star 4 used by these authors. The difference, 
in our data, in the R c magnitudes of the OT, resulting from the use of stars 1, 2 and 4, or only 
star 4 as a reference star, has an r.m.s of 0.027 magnitude. For data sets with a sufficiently large 
number of points and covering a time interval of at least 50 days on the power law decay part 
of the light curve (i.e., those of Sokolov et al. 1997, Pedersen et al. 1997, and our data) we 
made a second-order correction by a magnitude calibration offset, determined from the weighted 
average offset with respect to the power law decay part of the light curve. The corrections are 
0.027 ± 0.025, -0.020 ± 0.022 and 0.082 ± 0.031, for, respectively, our data (Tab. |), Sokolov et 
al. (1997) and Pedersen et al. (1997). We have assumed a conservative minimum error in the 
differential magnitudes of 0.05 magnitudes. 

The R c band light curve shows a maximum between May 10.2 and 10.8 UT, 1.3-1.9 days 
after the 7-ray burst. We fitted the rising and decaying parts of the light curve with power 
laws, Fr oc t a . For the fit to the rising part we use both the non-differential photometry and 
the differential data (without the second-order correction), assuming a minimum error of 0.10 
magnitudes. For the fit to the decaying part we only use the second-order corrected data. The 
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R c band light curve rises for —0.1 < log(t) < 0.20 with exponent a = 3.0 ± 0.4 (x 2 = 5.0/6). 
After the maximum (log(t) > 0.27) we find a = -1.141 ± 0.014 (% 2 = 49/20). Including in the fit 
the non-differential data and the other differential data, i.e., without the second-order correction 
(assuming a minimum error of 0.10 magnitudes) gives a = —1.132 ± 0.015 (x^ = 74.6/38). Since 
a hint of flattening of the light curve was reported (Pedersen et al. 1997), we have also fitted 
a power law plus a constant, Fr = Fn Q t a + C, to the decay part of the second-order corrected 
differential R c band light curve. We find a = —1.191 ± 0.024; the constant C corresponds to 
R c = 26.09 ± 0.36 (Xx = 43/19; note that this corresponds to R c = 26.30 ± 0.36 in the calibration 
of Pedersen et al. 1997). 

We also determined the optical fluence S of the OT (3000 A- 10000 A) taking for its spectral 
slope = -1 or -1/3 (see Sect. §. We find S = (1.20 ± 0.07) x 10~ 7 erg cm" 2 ((3 = -1), 
or S = (1.90 ± 0.12) x 10~ 7 erg cm" 2 ((3 = -1/3) for the fluence collected until day 110 after 
the event. Extrapolating the power law decay beyond day 110 we find a total optical fluence 
Soo = (1.36±0.08) x 10~ 7 erg cm" 2 ((3 = -1) and = (2.16 ±0.14) x 10~ 7 erg cm" 2 ((3 = -1/3). 

4. Broad-band optical flux distribution 

Between May 9.905 UT and May 10.033 UT the R c band flux increased by 0.137 ± 0.014 
magnitude per hour. In estimating the spectral energy distribution we corrected the May 9.9 UT 
observations for this brightening, i.e., we assumed that in all passbands the rate of brightening 
is the same, and reduced the magnitudes to a single epoch (May 9.93 UT). The May 10.98 UT 
observations occured at maximum light (see Fig. |TJ), so we have not applied any corrections. 
Corrections to the May 12 data were calculated using a power law fit to the decay part of the R c 
band light curve (see Sect. ||). The corrections are minor (< 0.04 mag) and marginally (< 0.02) 
affect the spectral slope. 

We fitted our B, V, R c , I c photometry with a power law, F v oc v", using the photometric 
calibration of Bessell (1979; table IV) (We chose to exlude the U band data to be able to compare 
our results directly with Sokolov et al., 1997. This exclusion does not affect the spectral slope 
significantly). Before the fits were made, the data were corrected for interstellar extinction Ay. 
Using the IRAS database, we have determined the 100 /im cirrus flux in the direction of the 
OT, Jioo < 0.08 MJy sr _1 (2a), corresponding to Ay < 0.01 (Laureijs et al. 1989). However, 
to compare our results with those of Djorgovski et al. (1997) we also used their value for the 
interstellar extinction (Ay = 0.08, corresponding to their Ab = 0.11). The effect of varying the 
interstellar extinction on the spectral slope is small; the maximum correction to the exponent 
(3 (with respect to Ay = 0.08) is -0.10 (for Ay = 0). We have also taken the B, V, R c and I c 
measurements of Sokolov et al. (1997) and determined (3, correcting for an extinction of Ay = 
0.08. In Figure |2] we have plotted the spectral slopes as a function of time, including the results of 
Djorgovski et al. (1997) and Metzger et al. (1997b). 
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5. Discussion 

The differential R c light curve shows that deviations from a smooth power law decay are 
moderate (r.m.s. = 0.15 mag; lower panel in Fig. [I]). The magnitude R c = 26.09 ± 0.36 for a 
possible underlying galaxy is consistent with the limits derived by Sahu et al. (1997b) and Pian 
et al. (1997), i.e., R c > 25.5. We note, however, that the introduction of a constant to the power 
law fit has not improved the x 2 substantially; we conclude that the evidence for a flattening of the 
light curve at late times is not yet convincing. 

The optical fluence Soo (3000A-10000A, = -1; Sect. |), is about 5% of the BATSE 20-1000 
keV and 20% of the WFC 2-26 keV GRB fluences. A similar comparison with the X-ray afterglow 
(Piro et al. 1997a; Sx = 7.3xl0 -7 erg cm -2 2-10 keV, from an extrapolation of the X-ray light 
curve for t > 27 s) shows that the optical afterglow fluence is ~ 20 % of the X-ray afterglow 
fluence. 

The optical spectrum of the OT becomes redder before reaching maximum light (Fig. 
Q); afterwards the spectral slope remains constant. We distinguish three phases of the light 
curve: phase I, the rising phase (a = 3.0 ± 0.4, (3 = —0.33 ± 0.17), phase II at maximum light 
(a ~ 0,/3 = -0.9±0.10), and phase III, the decaying phase (a = -1.141±0.014, j3 = -1.11±0.06). 
Wijers et al. (1997) give a relation between the spectral slope f3 and the temporal slope a for the 
simplest fireball remnant model (a forward blast wave with single Lorentz factor T; their Eq. 2) 
and for a 'beamed' fireball, i.e. T and the energy per unit solid angle E are functions of angle 
(their Eq. 4). For phase III we find very good agreement with a = —1.141 ± 0.014 implying 
f3 = —1.094 ± 0.009 for the 'beamed' case (marginal agreement with [3 = —0.761 ± 0.009 for the 
simple case). Phases I and II are not in agreement with the simple and beamed fireball models, 
that predict less steep or even inverted spectral slopes (3. However, these early phases might 
be explained by a continuous distribution of Lorentz factors (Rees & Meszaros 1997). Lower V 
(redder) material takes longer to sweep up a significant amount of external matter and hence its 
emission is delayed with respect to the higher V (bluer) material. The lower V material catches up 
gradually with the decelerated higher V material. Hence the rise of the light curve and evolution 
from blue to red can be explained with building up the emission from material with lower T, while 
the power law decay (phase III) sets in when most of the material has been mixed to a single T. 
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R-band lightcurve of GRB970508 
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Fig. 1. — Upper panel: the R c band light curve of GRB 970508. For the differential lightcurve 
(closed symbols) data have been used from Tab. (*), Sokolov et al. (1997; •), Pedersen et al. 
(1997; ■), Garcia et al. (1997; ▲), Schaefer et al. (1997; ♦) and Chevalier et al. (1997; x). The 
open symbol (o) represents non-differential data from Djorgovski et al. (1997), Castro-Tirado et al. 
(1998), Mignoli et al. (1997), Kelemen (1997), Fruchter et al. (1997) and Metzger et al. (1997c). 
The four open symbols (□) represent wide R filter data by Pedersen et al. (1997). The power law 
fit to the second-order corrected data with exponent a = —1.141 ±0.014 is represented by the solid 
line. Indicated by the dashed line is a fit to the second-order corrected data of a power law plus 
a constant, Fr = Fn t a + C, with exponent a = —1.191 ± 0.024. The constant C corresponds to 
a magnitude R c = 26.09 ± 0.36 Lower panel: the fit to the power law decay subtracted from the 
second-order corrected data of Tab. |l| (*), Sokolov et al. (1997; •) and Pedersen et al. (1997; ■). 
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Fig. 2. — The spectral slope f3 as a function of time in days after the burst (we used Ay = 0.08, 
see Sect. |j). Data are from: this paper (*), Sokolov et al. (1997; x), Djorgovski et al. (1997, •) 
and Metzger et al. (1997b, o). The maximum of the light curve occured between day 1.3 and 1.9 



Table 1: Log of observations of GRB 970508 
Date (UT) Telescope U B 



May 9.20 
May 9.93 
May 10.03 
May 10.98 
May 11.01 
May 12.03 
May 22.97 
July 4.19 
Aug 26.9 



WIYN 
WHT PF 
WHT PF 
WHT PF 
WHT PF 
WHT PF 
WHT AUX 
WHT AUX 
WHT PF 



20.61±0.25 
19.75±0.25 
19.97±0.25 



20.70±0.08 
21.21±0.08 



V 

21.05±0.05 
20.19±0.05 
20.71±0.05 



Rc 



21.25±0.05 
20.88±0.05 
20.46±0.05 
19.92±0.05 
19.77±0.07 
20.30±0.07 
22.04±0.07 
23.95±0.20 
24.57±0.07 



Ic 

20.49±0.18 
19.61±0.30 
19.89±0.18 
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Fig. 3.— WHT prime focus 3'.5 x 3'.5 R c band image of the field of GRB 970508. Indicated are 
Bond's variable (OT) and the 6 stars used for calibration. Note that Sokolov's star 3 is not indicated 
since we did not use it as a secondary photometric standard. The U band magnitudes of the four 
U band comparison stars, respectively, stars 2, 5, 6 and 7 are 19.92±0.25, 19.26±0.25, 18.99±0.25, 
and 17.61±0.25. 



